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Abstract 

An extended version of the non linear Walecka model, with p mesons 
and eletromagnetic field is used to investigate the possibility of phase 
transitions in hot (warm) nuclear matter, giving rise to droplet forma- 
tion. Surface properties of asymmetric nuclear matter as the droplet 
surface energy and its thickness are also examined. 
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1 Introduction 



In heavy-ion collisions, part of the hot nuclear matter produced can be de- 
scribed in terms of hadrons. The formation of highly excited composed nuclei 
in equilibrium with a gas of evaporated particles can be interpreted in the 
framework of hydrodynamics as two coexisting phases of nuclear matter, a liq- 
uid and a gas phase Q. During these reactions, phase transitions may occur 
depending on the temperature and densities involved The investigation 

of asymmetric nuclear matter is also of particular interest for problems in as- 
trophysics. In fact, neutron-star matter at densities between 0.03 fm -3 and 
nuclear matter density (0.17 fm -3 ) consists of neutron-rich nuclei immersed in 
a gas of neutrons ||. The size of the nuclei is determined by a competition 
between the surface energy and the Coulomb interaction. The use of ther- 
modynamical concepts in the study of possible phase transitions in the above 
problems is done with the underlying assumption that the time required for 
thermalization and chemical equilibrium is small. 

Hot nuclei, liquid-gas phase transitions and droplet formation in nuclear 
reactions as well as the surface properties of nuclear matter have already been 
extensively discussed in the literature in the context of non-relativistic models, 
namely within the framework of the Hartree-Fock (HF), Thomas-Fermi (TF) 
and Extended-Thomas-Fermi approximations (ETF) 0,11- fllPf . In partic- 
ular, in it is shown that the semi-classical TF approximation scheme is 
reasonably accurate at any temperature. 

Within the framework of relativistic models phase transitions in hot (warm) 
nuclear matter have also been investigated in infinite matter by imposing con- 
stant mesonic fields ||, at zero temperature for symmetric semi-infinite nuclear 



matter jLl]], at finite temperature for symmetric matter in the linear Walecka 
model and also in its non-linear form for symmetric and asymmetric 

matter [I3|. Surface properties of asymmetric semi-infinite nuclear matter 
have been investigated at zero temperature in |14| and in terms of a semi- 
classical treatment. In most of the above mentioned papers in which temper- 
ature effects have been taken into account, the finite temperature version of 
the liquid drop model is used in the investigation of the surface properties of 
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the arising droplet. 

In the present work, we study the possibility of droplet formation in a 
vapour system at finite temperature in the framework of the relativistic Walecka 
model with non linear terms (NLWM), which is known to describe adequately 
the bulk properties of nuclear matter. We include the Coulomb interaction 
and work in the Thomas-Fermi approximation. We determine the conditions 
for phase coexistence in a mult i- component system by building the binodals 
for several temperatures and two parametrizations of the non-linear Walecka 
model. These values determine the initial conditions which are used in solving 
numerically the coupled equations of motion obtained in the Thomas- Fermi ap- 
proximation at finite temperature. This semi-classical approximation contains 
the following quantal ingredients: the production of anti-particles is included 
and the distribution function of particles and anti-particles takes into account 
the Pauli Principle. 

In section 2 we obtain the equations of motion in the static case. The ther- 
modinamical potential in the framework of the Thomas-Fermi approximation 
is calculated in section 3. In section 4 the two-phase coexistence is discussed 
and in section 5 we present the numerical results. Finally, in the last section 
some conclusions are drawn. 



2 Extended Non Linear Walecka Model 

In what follows we describe the equation of state of non-symmetric matter 
within the framework of the relativistic non-linear Walecka model 111611, [117 



with the inclusion of p-mesons and the eletromagnetic field. The self interac- 
tion terms of the scalar meson were shown to be necessary in order to ade- 
quately describe nuclear properties []i~6| . Both the p meson and photons are 
incorporated to account, respectively, for the neutron excess in heavy nuclei 
and the eletromagnetic interaction between the protons JTTJ . In our lagrangian 



the 7r-meson field amplitude is not considered since we are not interested in 
pion-condensed states and hence, under the approximation we use here, all 
pion contributions vanish. 

In this model the nucleons are coupled to neutral scalar 0, isoscalar-vector 
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V^, isovector- vector meson fields and the electromagnetic field A^. The 
lagrangian density reads: 



C = ip 



where 



and 



% v = dfjVv - d v V„ (2) 
= df+K - d u b p - g p {b^ x b„), (3) 

Ffj,v = d^A u - d v Afj,. (4) 



The model comprises the following parameters: three coupling constants g s , 
g v and g p of the mesons to the nucleons, the nucleon mass M, the masses of 
the mesons m s , m v , m p , the eletromagnetic coupling constant e = y^§f and 
the self-interacting coupling constants k, A and £. We have used two sets of 
constants, respectively identified as MS taken from |2j and as NL1 taken from 
22] , with Cf = gfM 2 /mf, i = s, v, p, and they are displayed in table 1. 



All meson masses must be specified along with the coupling constants. It is 
true that only the rations of couplings to the masses are necessary in infinite 
matter in the linear Walecka model. When non - linear terms are included, 
the statement looses its validity. 

From the Euler-Lagrange formalism, we obtain the coupled equations of 
motion for the scalar, isoscalar-vector, isovector-vector, eletromagnetic and 
nucleon fields, respectively given by: 

{dl - V 2 + m 2 )0 = g sPs - ^ - (5) 
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{d 2 _ v 2 + ml)V» = g v f - |<rf(V») 3 , 



(6) 



{d 2 t - V 2 + m 2 ,)^ = + x + x P*), (7) 



(5* - V 2 )^ = 



and 



idtijj = 



a • HV - a,V - |r 3 b - e2^A) 



(9) 



where the scalar density p s and the baryonic current densities are defined as 

p s =< 4>tp >, 

f =< jrf$ >, 
f =< ^f^J >, 

fern =< Wt^ 1 )* > 

and 63 = (6°,b). In the static case there are no currents in the nucleus and 
the spatial vector components V, b and A are zero. Therefore, the equations 
of motion become: 



V 2 = m 2 + ^0 2 + ^A0 3 - g sPs , 



(10) 



V 2 V = m 2 v V + y g 4 X-9vPB, 



V\ = m 2 p b - ^p 3 , 



(11) 
(12) 



V 2 A = -ep p , 



(13) 



where P b = P P + Pn and p 3 = p p — p n are the barionic densities, and Pp and p n 
are the proton and neutron densities. 



5 



3 The Thomas— Fermi Approximation 



The present work is based on the semi-classical Thomas-Fermi approxima- 
tion. In this approach the energy of the nuclear system with particles and 
anti-particles, described respectively by the one-body phase-space distribution 
functions n + (r, p, t) and n_(r, p, t), at position r time t and with momentum 
P 

, ,x f n p ±(r,p,t) \ 

n±(r,p,t) = , (14) 

\ n„ ± (r,p,t) J 

is (only the nuclear matter contribution and interaction terms) 

En = 7Tr J (n+(r, p, t)e+ + n_(r, p, t)e_) (15) 

where 

- = | Cp± ° ) , e l± = v / (PTV 4 ) 2 + (M-^0) 2 ±V iO , / = /;■//. 
e n± 



with 



V p o = ^«Vb + y&o + eA , V n0 = g v V - y6 , 



V p = g v V + g p b + eA , V„ = #„V - , 

are the classical effective one-body Hamiltonian for particles (+) and anti- 
particles (-) since particles and anti-particles have opposite baryonic charge and 
7 = 2 refers to the spin multiplicity. We can also work with the distribution 
function for particles at position r, instant t with momentum p, /+(r, p, t) = 
n + (r,p,t) and the distribution function for antiparticles at position r, instant 
t with momentum — p, /_(r, p, t) = n_(r, —p,t) so that 

En = 7Tr J ±0± (f + (r,p,t) h + -/_(r,p,i) h.) 

where 

h± = ±e±(r, ±p,t) 

±v/(p-V p ) 2 + (M-^0)2 + V p o 

± V / (p-V„) 2 + (M-^0)2 + V nO 
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The classical entropy of a Fermi gas is given by 

+ 1.1(1 - /,+) + (/,+ « /,_; 



i=p,n 



(27T) 



(16) 



and the thermodynamic potencial is defined as 



Q = E - TS - Vi B i 



(17) 



i=p,n 



where B„ , B n are, respectively, the proton and the neutron number: 



B, 



d 3 rp;(r,t), A = 7 



d 3 p 
(2^ 



(f i+ -fi-), i = p,n, (18) 



fii is the chemical potential for particles of type i and T is the temperature. 
For a system in equilibrium, the distribution functions should be chosen to 
make the thermodynamic potencial Q stationary and hence 

1 



fi±{r,p,t) 



(19) 



1 + exp[(e T ^)/T] ' 

where Vi = fa — Vjo is the effective chemical potential, e = ^p 2 + M* and 
M* = M — g s (f) is the effective nucleon mass. In the static approximation 
V± = 0. 

From the above expressions we get for (|T7|) 

d 3 rd 3 p 



tt = 7 Tr 



(2tt) j 



(f+(r,P,t)h + - /_(r,p,t)fr_ 



+i | d 3 r [(V0) 2 - (W ) 2 - (VM 2 - (VA ) 2 ] 



+- / d 3 r 



m 2 2 + + §A0 4 - - %tgX - ™y Q 



3! 



4! 



- 7 T^/d 



3 d'p 

r 

(2tt) 3 



- + ln(l + e -(-^) + ln(l + e^ + ^ T ) 

(20) 
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Equation ( PU[ ) can then be written in the form 

O = j d 3 r Q [(V0) 2 - (W ) 2 - (V6 ) 2 - (VA ) 2 ] - V ef ^j 



with 



mW + §/# 3 + |A0 4 - m% 2 - Ifc'ltf ~ mjftg 



4! 



E / 7^3 N 1 + e" (£ -^ )/T ) + ln(l + e~^)/ T 
The fields that minimize Q satisfy the equations 



dV f 



- m s4> ~ 2^ - ^j A< ^ + 9sP. 



dV f 



ef 



dV ef _ m 2 h 9 Pn 



9V F 



e/ 



Mr 



-ep P , 



(21) 

(22) 
(23) 
(24) 
(25) 



where Pb, P3 were defined at the end of section 2, p p has been defined in (jl 
and 

d 3 p M* 



= 7 



i=p,n 



(2vr) 3 e 



Comparing eqs. (|TU|-|T^) with eqs. 



we see that 



d 2 (f> 2d0 
dr 2 r dr 



90 ' 



V 2 6 



V 2 A 



dV , 2dV dV ef 
H 



dr 2 



+ 



r dr 
2d6 



(26) 
(27) 
(28) 

dr z r dr oAq 
These coupled differential equations are solved numerically and all relevant 

quantities (e.g. effective mass, densities, pressure), which depend on the fields 

are calculated. 



dr 2 r dr 
d 2 A 2 dA 



dV ' 

dVef 
db ' 

dVef 
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4 Two— phase coexistence 



In order to obtain the initial conditions for the program which integrates the 



differential equations (|26| ~ p9|) we determine the conditions under which two 
distinct phases can coexist in infinite matter. In the mean field approximation 



the meson fields are replaced by their expectation values pi 18 



=< <p >= 00, 



(30) 



V° =< V° >= V , 



(31) 



b u =< W >= b , 



(32) 



and the eletromagnetic field vanishes. The substitution of the above expres- 
sions in eqs. fllPf) , flllD and ( |12|) yields 



2m 2 /" 6 



2 ^ ,3 , 9s 



mt 



(33) 



V a = - 



M 

6ml 



T 7"3 I 9V 

V + —PB 

mi 



(34) 



7 9 P 



(35) 



The thermodynamic quantities of interest are given in terms of the above 
meson fields. They are the energy density: 

£ = 2^ E / P 2 d P \/p 2 + M* 2 (f l+ + /,_) 

i=p,n 



+ 2 V ° + 8 Vo+ 8mf 3 + 2 00 + 6 00 + 24 00 ' 



(36) 



the pressure: 



7 

67T 2 



i=p,n 



p 4 dp 
\Jp 2 + M* 



(/*+ + /i- 
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the entropy density: 



<S = + P - UpPp- finPn), (38) 



and the proton fraction: 



Y p = (39) 
Pb 



A thorough study of the possibility of phase transitions in hot, asymmetric 
nuclear matter is done in p(J and In figure 1 we plot the pressure in 



terms of the baryonic density ps for each proton fraction and for T = 10 MeV 
obtained with the MS constants. This figure is slightly different from figure 4 
of @, which is reproduced when m s = 500 MeV. Similar behaviours are found 
for T = 5 MeV and also with the NL1 constants. 

We have made use of the geometrical construction |2(J in order to obtain the 
chemical potentials in the two coexisting phases for each pressure of interest. 
As an example, we show p p and p n in function of the proton fractions in figure 
2 again for T = 10 MeV and MS constants. 

In a binary system 

dp p \ ^ q d ( ^^ n \ < o (40) 

OYp / rp p \ 9Yp J rp p 

known as diffusive stability, which reflects the fact that in a stable system, 
energy is required to change the proton concentration while pressure and tem- 
perature are kept constant. In order to obtain the binodal section which con- 
tains points under the same pressure for different proton fractions, we have 
used the conditions above and simultaneously solved the following equations: 

P = P(v p ,v n ,M*), (41) 



P = P(u' p ,u' n ,M*'), (42) 



t n{v p ,v n ,M*)= t jL i (i/ p ,i/ n ,M*') ,i=p,n (43) 
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^0O+^30o+^O = P S (^^n,M*) (44) 



and 

^0o' + + = ^, M*'). (45) 



The binodal sections for the MS and the NL1 constants and temperatures 
equal to 5 and 10 MeV are plotted, respectively, in figure 3 and in figure 4. For 
certain values of proton and neutron chemical potentials, the system may be 
at the same pressure with different densities and proton concentrations, which 
allows for the possibility of phase transitions. For the sake of completeness, 
we also show in tables 2, 4 and 6 some of the points taken from the binodal 
sections. The results we have chosen as input to the code which solves the 
differential equations (^6| - ^9|) are displayed in the last three columns of these 
tables. 



5 Numerical Results 



Solving numerically the set of coupled equations (p6|-[29|) is not trivial. The 
main problem is related with the boundary conditions which have to be set 
within the droplet. To understand better this statement, please refer to the 
Appendix, where the formulae are simplified to the T = case and this prob- 
lem becomes clear. At first, we have tried to use the code COLSYS ||21|| , as 



suggested in [f2^] and [0, but we have obtained satisfactory results only for 



symmetric nuclear matter | 23| . We have finally opted for another code, written 
with the help of the Gears stiff integration method, which uses as input the 
temperature, the size of the mesh, boundary conditions and initial conditions. 
The chemical potentials are output and can be fitted in accordance with the 
size of the mesh. In the Appendix the boundary and initial conditions we have 
used are explicitly written. 

The radius R ma x (see Appendix) fixes the neutron chemical potential calcu- 
lated at the last mesh point. The proton chemical potential is fixed by another 
(inner) mesh point, which depends on the difference desired between both po- 
tentials. If fi p = fi n (symmetric nuclear matter), R ma x fixes both chemical 
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potentials at the last mesh point. We have considered that convergence has 
been achieved when the baryonic density does not vary more than 0.5 per cent 
and the chemical potentials are close to the ones obtained from the binodal 
section for a certain proton fraction. 

As an example, in figure 5 we plot the fields which are solutions of the 
coupled equations for the MS parametrization at T=10 MeV with the initial 
conditions given in table 2 for Y p = 0.3. The corresponding barionic density 
is plotted in figure 6 and represents a droplet of the liquid phase (small r) 
in the background of the vapour phase (large r). Similar density profiles are 
obtained for all other proton fractions, except for Y p = 0.5, when the curves 
for proton and neutron densities coincide. The region at the surface with extra 
neutrons is known as neutron skin. Some quantities of interest to study the 
surface properties are the two squared-off radii R n and R p , defined in || as 

^Ftmax 

Pn{r)dr = p n ,iR n + Pnj(Rmax ~ Rn) , (46) 



and 



I 

Jo 



p p {r)dr = p P;i R p + p p j(R m ax ~ Rp), (47) 



where pi refers to the liquid density, p/ to the gas density, R ma x is the size 
of the mesh for which convergence is achieved, and the neutron skin thickness 
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= R n - Rp. (48) 

These quantities are computed for the droplet solutions we obtain and given 
in tables 3, 5, 7. 

The droplet surface energy and thickness are obtained from the free energy 
of a system with a fixed number of particles B = B p + B n , in which a droplet 
of arbitrary size grows in the background of the vapour phase. Within the 
small surface thickness approximation it reads ||12j|: 



Anr 2 dr 



2 _ /dA x 2 

\dr J \ dr J \ dr J \ dr 



PpB p ~\~ p n B, n 
(49) 
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For droplets with radius R and small surface thickness, 



F(R) =47caR 2 -CV + n P B p + ^ n B n , (50) 

where C is a constant and V is the volume of the system. The surface energy 
per unit area of these droplets is 



DC 



a = dr 



d<A 2 fdV \ 2 fdb \ 2 fdAt 



dr \ dr \ dr \ dr 



(51) 



The surface thickness t is defined as the width of the region where the density 
drops from 0.9pBo to O.lpsO) where pso is the baryonic density at r = 0. 
According to fi~6fl , for T = 0, o should be of the order of 1.25 MeV fm -2 and 



t of the order of 2.2 fm. 

In table 2 some points taken from the binodal section for the MS constants 
at T = 10 MeV are explicitly written. In table 3 results found for the proton 
fractions at r = (Y p (i)), chemical potentials, the surface energy, its thickness, 
the size of the mesh for which convergence is achieved, R ma x, Rp and the 
neutron skin thickness are displayed. Notice that there is a small discrepancy 
between the proton and neutron chemical potentials given in table 2 and the 
ones displayed in table 3. This is due to finite size effects and the inclusion of 
the Coulomb interaction. In tables 4 and 5 again the points obtained from the 
binodal section at T = 5 MeV and the respective droplet solutions are shown 
for the MS set of constants while in tables 6 and 7 the NL1 constants are used. 

At this point some comments are in order. The size of the mesh (R max ) 
given in tables 3, 5 and 7 are the smaller values for which there is convergence 
for a given density at r = and a given proton fraction Y p within the accuracy 
of the present numerical calculations, i.e., ±0.5%. We have chosen to compare 
data corresponding to the same value of the proton fraction at r = because 
this parameter is independent of the propeties of the surface. A larger mesh 
size would converge to a larger droplet with the same values for the proton 
fraction Y p and the density at r = 0, and the same chemical potentials. For 
a detailed explanation on the introduction of the Coulomb field, please refer 
to the Appendix. Before drawing our conclusions, we would like to emphasize 
that in our calculations, the proton and neutron numbers are never fixed. They 
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are just consequence of the results for the fields and densities obtained from 
the convergent solutions of the differential equations. 



6 Conclusions 

We first examine the behaviour of the total baryonic density. From figure 6, 
one can see that it falls from the initial liquid density to the vapour density, 
which is very small, but different from zero, as expected. One can compare 
this figure with the densities presented in table 2 for infinite nuclear matter 
with Y p = 0.3. Concerning the neutron and proton densities, figure 6 shows 
the same profile obtained in |TD[ and H with non-relativistic models. 



The proton fraction in the vapour phase is smaller than in the liquid phase 
(Y p (i)), except for the symmetric nuclear matter, when it remains unaltered. 
This can be seen from tables 2, 4 and 6 and confirmed for the droplet solutions. 
This fact can be interpreted as a nucleus with a given proton concentration 
(the phase of higher density) in equilibrium with a gas of drip nucleons, mostly 
neutrons (the lower density phase) with a much smaller proton concentration. 

From tables 3, 5 and 7, one can check that the surface energy a increases 
with the initial proton fraction and its thickness t decreases. In fact, the 
larger the proton fraction the less important is the contribution from the b 
field in the a calculation as can be seen from fl5lf) . For symmetric nuclear 



matter, the results are compatible with the ones suggested in |L6) and for 
asymmetric matter the surface thickness results are comparable with the ones 
presented in \TE\. One can also see that the larger the proton fraction, the 



smaller the size of the mesh for which convergence is achieved. This may be 
due to the decrease of neutron-proton asymmetry and therefore, the increase 
of the droplet binding. Some conclusions with respect to the temperature 
dependence of the droplet solutions can be drawn comparing tables 3 and 5. 
First, it is easy to note that the size of the mesh must be larger for higher 



temperatures. A similar statement was made in ||13|| . Also, the surface energy 
decreases with the increase of temperature, while the surface thickness is larger 
for higher temperatures. The decrease of a with the temperature is easily 
understood from fl5T|), because for higher temperatures the fields decrease more 
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smoothly and spread out over a larger distance at the surface and, therefore, 
their derivatives are smaller. 

In the same tables the squared-off proton radius and the neutron skin 
thickness are also shown. Both quantities decrease with the increase of the 
proton fraction at r = 0. The numbers we have obtained can be compared 
with the ones found in M. These behaviours could be a consequence of the 
increase of the droplet binding. The neutron skin thickness is a quantity which 
is larger for lower temperatures, except for symmetric nuclear matter, when it 
is zero, independently of the temperature considered. 

Concerning the importance of the Coulomb interaction and its consequences 
in the droplet formation, one can see, from figure 6, that the proton and neu- 
tron densities are indeed modified by the eletromagnetic field, as pointed out 
in Jl| . Nevertheless, for the same value of Y p at r = 0, the profile of the sur- 
face with eletromagnetic field is almost unaltered with respect to the results 
obtained without the Coulomb interaction. This fact is reflected in the results 
displayed in table 3. a, which were obtained with the same input parameters 
as the ones shown in the third line of table 3, but without the inclusion of the 
eletromagnetic interaction. The surface energy a, the proton radius R p and 
the neutron skin are slightly different when they are calculated with and 
without the Coulomb potential while t is practically the same. Two effects of 
the eletromagnetic field in the present calculation is to decrease the number 
of particles in the droplet, since the central density becomes smaller, and to 
increase the proton radius Rp. 

The conclusions drawn above are independent of the sets of parameters 
used. 

To the best of our knowledge, this is the first time that the eletromagnetic 
interaction is taken into account within the framework of a relativistic model 
in order to calculate surface properties. We believe that a more systematic 
study of the importance of the Coulomb field for various temperatures and 
proton fractions has still to be made and this problem will be tackled in a 
forthcoming work. 
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Appendix 



For T = the distribution functions are simply given by: 



and, hence, 



1 

~2 

"7 



hi — hi + , 



E 



i=p,n 



3! 

d 3 p 
(2^ 



4! " u 4! 

hifi ~\~ PpPp ~\~ Pnprii 



Ps = l 



E 



d 3 p M* 



fi 



7 ^ f kMr) _ 2 . M 



E 



(2vr) 3 e ' 2tt — ,„ 

i=p,n 4 ' i=p,n u 



p dp- 

e 



., (27r) 3,; 6tt 2 

Minimization of f2 with respect to k F i(r), i = p,n, gives rise to the following 
conditions 

4 P (r) - yfk 2 Fp (r) + M* 2 (r) - g v V - |& - eA,) = 



and 



k 2 Fn (r) Ln - ^Jk Fn (r) + M* 2 (r) - g v V + |& ) = 0. 



We obtain k Fp {r) = and k F n(f) = or, for k Fp {r) and k Fn (r) different from 
zero, 



= J k Fp (r) + M ( r ) + 9vV + g p b + eA , 



(52) 



Pn = \/kl n {r) + M* 2 (r) + ^V Q - g p b . (53) 

The value of k Fp {r) and k Fn {r) is obtained inverting these last two equations. 
The discontinuity on the values of k Fp (r) and k Fn (r) has to be taken into 
account in the code that solves the differential equations (|22"|-f25|) . Outside the 
droplet k Fp (r) and k Fn (r) are zero and the mesons are free. 
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In our code, the boundary conditions are given by 

d0, s dV db , s dA , , 

— (r = 0) = — — (r = 0) = -^(r = 0) = — — (r = 0) = 
dr dr dr dr 

and for r = R ma x, where R max is the size of the mesh, 

# / 1 W n 

u -' - 1 4nax 

with similar equations for Vq an d &o an d 

dr Rmax 

or, considering the electron screening, v4 is zero at the last point of the mesh, 

Both boundary conditions give similar results. 

As initial guesses for the meson fields we have used Fermi like functions 
such as 

Mr) = h , 

1 + exp(m s (r - 0.8R max )) ' 

for the scalar field. The values for r = 0, i.e. <f> , V o and b were obtained 

from the binodal section. An initial guess for the eletromagnetic field is the 

field of a homogeneous spherical distribution of protons which at r = is of 

the order of 1 x 10" 2 in units of nucleon mass. The convergence is obtained 

for different initial guesses except if the initial guesses are much stronger than 

the value refered above. We also suppose that the droplets are formed in an 

electrically neutral enviroment, as we would find in neutron stars. 
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Figure Captions 



Figure 1. The pressure in terms of the baryonic density is plotted for 
each proton fraction for the MS set of constants. The first curve on the 
left relates to Y p = 0, the second one to Y p = 0.1 the last one to 
Y p = 0.5. The temperature is T = 10 MeV. The pressure is given in 
MeV/fm 3 and the density in fm~ 3 . 

Figure 2. The proton (lower curve) and the neutron (upper curve) 
chemical potentials are shown in function of the proton fraction for the 
pressure of 0.12 MeV/fm 3 . Again this graph is plotted for MS constants 
and T = 10 MeV. 

Figure 3. Binodal section for T = 10 (dashed line) and T = 5 MeV 
(solid line) with the MS constants. The pressure is given in MeV/fm 3 . 

Figure 4. Binodal section for T = 10 (dashed line) and T = 5 MeV 
(solid line) with the NL1 constants. The pressure is given in MeV/fm 3 . 

Figure 5. <P,Vo,Aq and b (in this order, from top to bottom) are shown 
in terms of r. The fields are given in nucleon mass units and are obtained 
with the input values of table 2 for Y p = 0.3. 

Figure 6. From top to bottom the density profiles for the baryons ps(r), 
the neutrons p n (r) and the protons p p (r) in fm~ 3 are plotted for the same 
case as in Fig. 5. The solid curves were obtained with the inclusion of 
the eletromagnetic field and the dashed ones without it. 

Table Captions 

Table 1. Sets of parameters used in this work. All masses are given in 
MeV. The * is a reminder that the authors of use a different scalar 
meson mass. 

Table 2. Results obtained from the binodal section built with the MS 
constants for T = 10 MeV. The pressures are given in MeV/fm 3 , the 
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densities in fm 3 , the chemical potentials in MeV and the fields in units 
of nucleon mass. 

Table 3. Output results given by the solution of the coupled differential 
equations with the MS constants and T = 10 MeV. Index i refers to 
r = 0. The effective chemical potentials are given in units of nucleon 
mass, the chemical potentials in MeV, the surface energy in MeV/fm 2 
and the surface thickness in fm. R max is the size of the mesh for which 
convergence is achieved and it is given in fm. R p and 6 are also given in 
fm. 

Table 3. a Output results given by the solution of the coupled differential 
equations with the MS constants and T = 10 MeV without the inclusion 
of the eletromagnetic field. The units are the same as in Table 3. 

Table 4. Results obtained from the binodal section built with the MS 
constants for T = 5 MeV. The units are the same as in Table 2. 

Table 5. Output results given by the solution of the coupled differential 
equations with the MS constants and T = 5 MeV. The units are the 
same as in Table 3. 

Table 6. Results obtained from the binodal section built with the NL1 
constants for T = 5 MeV. The units are the same as in Table 2. 

Table 7. Output results given by the solution of the coupled differential 
equations with the NL1 constants and T = 5 MeV. The units are the 
same as in Table 3. 
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Table 1 - parameters 



Force 


[ref] 




ci 


ci 


k/M x 10" 3 


A x 10- 3 




M 


m s 


m v 


m p 


MS 


[1] 


374.770 


260.570 


106.91 


3.0809 x 10 3 


8.106 x 10 3 


0.02364 


939. 


550.00* 


783.00 


770.00 


NL1 


[12] 


373.176 


245.458 


149.67 


2 gl 2.4578 


- 6 g\ 3.4334 


0.0 


938. 


492.25 


795.36 


763.00 



Table 2 





V 






Pb 


0o x 10" 2 


V 00 x lO" 2 


b x 10~ 2 


0.10 


0.38 


873.90 


941.36 




0.08 


2.01 


1.30 


-0.35 


0.18 










0.10 


2.69 


1.76 


-0.39 


0.20 


0.35 


874.66 


940.86 




0.11 


2.82 


1.85 


-0.38 


0.08 










0.07 


1.82 


1.17 


-0.33 


0.25 


0.27 


877.55 


939.04 




0.13 


3.15 


2.07 


-0.36 


0.05 










0.05 


1.35 


0.85 


-0.25 


0.30 


0.18 


882.16 


936.42 




0.14 


3.44 


2.28 


-0.33 


0.03 










0.03 


0.83 


0.51 


-0.16 


0.35 


0.10 


888.56 


932.55 




0.15 


3.69 


2.44 


-0.26 


0.03 










0.01 


0.41 


0.25 


-0.08 


0.40 


0.06 


896.36 


927.37 




0.16 


3.87 


2.57 


-0.19 


0.06 










0.01 


0.21 


0.12 


-0.04 


0.45 


0.03 


904.90 


920.77 




0.17 


3.98 


2.64 


-0.10 


0.19 








4. 


x 10~ 3 


0.11 


0.07 


-0.01 


0.50 


0.03 


913.20 


913.18 




0.17 


4.01 


2.67 


-l.OxlO* 4 


0.50 








3. 


x 10~ 3 


0.09 


0.05 


-0.1xlO~ 4 



Table 3 





u p (i) 




P P 




a 




Rmax 


R P 


e 


0.15 


0.75170 


0.78836 


876.18 


935.30 


0.22 


5.03 


16.70 


12.46 




0.75 


0.23 


0.70667 


0.73960 


879.62 


935.25 


0.35 


4.72 


15.73 


11.63 




0.71 


0.30 


0.66064 


0.68956 


884.29 


933.99 


0.58 


3.36 


10.49 


7.50 




0.59 


0.36 


0.62472 


0.64845 


892.46 


932.22 


0.70 


3.21 


9.44 


6.75 




0.47 


0.38 


0.61380 


0.63467 


896.94 


931.08 


0.76 


3.02 


8.39 


5.96 




0.38 


0.50 


0.61306 


0.61328 


913.20 


913.20 


1.21 


2.46 


7.34 


3.72 




0.0 


Table 3.a 




v p {i) 


v n {i) 


Pp 




a 




Rmax 


R p 


e 


0.30 


0.65717 


0.68534 


882.72 


933.99 


0.60 


3.36 


10.49 


7.39 


0.62 
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Table 4 





V 






Pb 


0o x 10" 2 


Vbo x lO" 2 


b x 10~ 2 


0.15 
0.05 


0.52 


868.12 


948.93 


0.12 
0.09 


3.01 
2.28 


1.97 
1.48 


-0.48 
-0.45 


0.20 
0.02 


0.37 


871.18 


946.76 


0.13 
0.07 


3.26 
1.92 


2.15 
1.23 


-0.46 
-0.39 


0.25 
6. x 10~ 3 


0.19 


875.63 


943.75 


0.14 
0.05 


3.49 
1.46 


2.31 
0.93 


-0.41 
-0.30 


0.30 
2. x 10" 4 


0.06 


881.47 


940.25 


0.15 
0.02 


3.73 
0.62 


2.47 
0.38 


-0.36 
-0.12 


0.35 
2. x 10" 4 


0.02 


888.80 


936.36 


0.17 

4. x 10" 3 


3.96 
0.13 


2.63 
0.08 


-0.29 
-0.03 


0.40 
2. x KT 3 


6. 

xl0~ 3 


897.55 


930.92 


0.18 
1. x 10~ 3 


4.12 
0.04 


2.74 
0.02 


-0.20 
-0.01 


0.45 
0.03 


I. 

xl0~ 3 


906.87 


924.01 


0.18 
3. x 10" 4 


4.22 
0.01 


2.81 
0.01 


-0.10 

-2. x 10~ 3 


0.50 
0.50 


5. 
xlO" 4 


915.82 


915.89 


0.18 
I. x 10~ 4 


4.25 
3. x 10~ 3 


2.83 
2. x 10~ 3 


0.0 
0.0 



Table 5 



Y p (i) 


u p (i) 




P P 


pn 


a 


t 


Rmax 


Rp 


e 


0.10 


0.75486 


0.79229 


873.60 


936.92 


0.22 


3.27 


12.59 


8.63 


1.16 


0.22 


0.70784 


0.73992 


879.47 


936.02 


0.42 


2.94 


10.49 


7.02 


0.98 


0.28 


0.67424 


0.70386 


882.57 


935.11 


0.58 


2.77 


9.23 


6.15 


0.81 


0.38 


0.91239 


0.63308 


899.08 


934.24 


0.92 


2.49 


7.34 


4.74 


0.50 


0.50 


0.59357 


0.59381 


915.80 


915.80 


1.35 


2.33 


5.24 


3.69 


0.0 
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Table 6 



Y P 


V 






Pb 


0o x 10" 2 


Vbo x lO" 2 


b x 10~ 2 


0.10 
0.11 


0.29 


885.36 


946.59 


0.07 
0.06 


1.90 
1.82 


1.11 
1.06 


-0.35 
-0.34 


0.15 
0.07 


0.27 


885.91 


946.13 


0.07 
0.05 


2.10 
1.61 


1.24 
0.92 


-0.35 
-0.32 


0.20 
0.06 


0.24 


886.59 


945.65 


0.09 
0.05 


2.43 
1.51 


1.45 
0.86 


-0.35 
-0.31 


0.25 
0.03 


0.18 


888.97 


944.02 


0.10 
0.04 


2.73 
1.23 


1.64 
0.69 


-0.33 
-0.26 


0.30 
0.01 


0.10 


892.61 


941.56 


0.11 
0.03 


3.02 
0.85 


1.83 
0.47 


-0.30 
-0.19 


0.35 
3. x 10~ 3 


0.04 


897.60 


938.35 


0.12 
0.01 


3.35 
0.36 


2.05 
0.19 


-0.25 
-0.08 


0.40 

5. x 10~ 3 


0.01 


904.32 


933.97 


0.13 

3. x 10~ 3 


3.67 
0.10 


2.26 
0.05 


-0.18 
-0.02 


0.45 
0.05 


4. 
xlO" 3 


912.06 


928.06 


0.14 
8. x 10" 4 


3.89 
0.03 


2.40 
0.01 


-0.10 

-5. x 10~ 3 


0.50 
0.50 


1.65 

xl0~ 3 


920.42 


920.46 


0.14 
3. x 10" 4 


3.97 
0.01 


2.45 
0.01 


-2. x 10" 4 
0.0 



Table 7 



Y P {i) 


v p [i) 






/in 


a 




Fornax 


R p 





0.15 


0.81150 


0.83765 


885.55 


937.21 


0.21 


3.55 


13.64 


9.56 


1.02 


0.24 


0.77241 


0.79539 


889.32 


936.84 


0.34 


3.27 


12.59 


8.84 


0.86 


0.32 


0.72522 


0.74488 


895.33 


936.27 


0.56 


3.00 


11.54 


8.26 


0.64 


0.38 


0.68604 


0.70168 


902.60 


934.54 


0.81 


2.65 


9.44 


6.80 


0.44 


0.42 


0.65833 


0.67021 


909.25 


932.60 


1.00 


2.51 


8.39 


6.15 


0.28 


0.45 


0.63661 


0.64454 


914.91 


929.43 


1.17 


2.50 


7.34 


5.45 


0.15 



25 



0.75 



0.5 



0.25 







-0.25 



-0.5 



-0.75 








0.03 



0.06 0.09 0.12 

Baryonic density 



0.15 0.18 



940 



920 



> 

CD 



CO 

-I—' 

£ 900 

Q_ 

"co 
o 

"E 

CD 
.£= 

o 



880 



860 








0.1 



0.2 



0.3 



0.4 



0.5 



Yp 



0.35 



0.3 - 



0.25 - 



0.2 - 



0.15 - 



0.05 - 




0.1 - 



0.04 



0.03 



0.02 



CD 



0.01 







-0.01 



0.15 



c/) 

CD 



0.12 



0.09 



to 

CD 
T3 



0.06 



0.03 







